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Ultrasonic transducers are a key element that governs the performances of 
both generating and receiving ultrasound in an ultrasonic measurement 
system. Electrical impedance is a parameter sensitive to the environment of 
the transducer; it contains information about the transducer but also on the 
medium in which it is immersed. Several practical applications exploit this 
property. For this study, the model is implemented with the VHDL-AMS 
behavioral language. The simulations approaches presented in this work are 
based on the electrical Redwood model and its parameters are deduced from 
the transducer electroacoustic characteristics. 
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1. INTRODUCTION 

Ultrasound systems are widely used. They find many applications in engineering, medicine, biology, 
and other areas. The one indispensible part in these systems is the transducer. These will use the properties of 
magnetostrictive or piezoelectric materials to convert electrical energy into ultrasonic mechanical energy [1]. 
Piezoelectric materials have the advantage over other systems of having good performance and being 
available in very diverse geometries. The electromechanical interaction of piezoelectric transducer, 
represented by electrical equivalent circuits, was first introduced by Mason [2]. He proposed an exact 
equivalent circuit that separated the piezoelectric material into an electrical port and two acoustical ports 
through the use of an ideal electromechanical transformer. The problems with the model are that it required a 
negative capacitance at the electrical port. Redwood [3] improved this electromechanical model by 
incorporating a transmission line, making possible to extract useful information on the temporal response of 
the piezoelectric component. 

The Electrical impedance is a parameter sensitive to the environment of the transducer. Several 
practical applications exploit this property[4, 5]. The measurement of the impedance makes it possible, for 
example, to detect the physical or structural modifications of the medium due in particular to damage. This 
approach is used for non-destructive testing to monitor the condition of structures such as aging and 
corrosion [6, 7]. Frequency analysis of the impedance makes it possible to precisely locate the resonance 
zone of the transducer. This location can be exploited to control and stabilize the operating frequency of high 
power systems such as ultrasonic welding devices [8]. Real-time knowledge of the electrical impedance also 
makes it possible to determine and optimize the power emitted by a transmitter or the sensitivity in reception 
of a piezoelectric sensor [9]. 
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Different approaches were proposed to predict the piezoelectric transducer behavior such as a 
numerical resolution of piezoelectric equations. Another approach is based on the equivalent electrical circuit 
simulation using an electric simulator like SPICE (Simulation Program with Integrated Circuit Emphasis) 
[10, 11]. However SPICE presents some limitations SPICE is on a continuous basis it cannot support discrete 
representations, and as a result, it is not suitable for mixed modeling [12, 13]. VHDL-AMS (Very High speed 
integrated circuit Hardware Description Language Analog and Mixed Signal) is a high-level language that 
enables numerical and analog simulations, while giving the possibility to simulate systems with different 
physical magnitudes: mechanical, thermal and electrical. The use of an analysis tool, such as the VHDL- 
AMS behavioral description language, can be a solution to the limitations caused by the use of the SPICE 
simulator. 


2. THEORY OF PIEZOELECTRIC ULTRASONIC TRANSDUCER 

The piezoelectric material is the principal element in an ultrasonic transducer. The piezoelectric 
materials have the advantage compared to the other systems of presenting good performances and to be 
available in very diverse geometries. These materials are generally appeared as a disc, a ring or plate. Our 
piezoelectric element has the characteristic to vibrate in thickness mode, on only one direction axis (z). This 
modelling study is thus limited to one geometrical dimension (1-D), with the boundary conditions at the 
acoustic ports Figure 1 [14]. 

F = -AT ( 1 ) 

Where A is the area of the transducer and T is the internal stress. 




Figure 1. Piezoelectric plate of thickness e and its representation as a three port system 


Let’s consider e the thickness of the piezoelectric plate and V 3 is the excitation voltage. Fi and F 2 
presented the forces transmitted to the propagation medium on the front and back face of the transducer, vi 
and V 2 are the acoustic particles velocities at the front and the back faces of the transducer. 

In these conditions, the transfer matrix form (2) which describes the global behavior between the 
electric excitation port and the two acoustic ports. 
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vi and V 2 (m/s) are the acoustic particle velocities at the front and the back faces of the plate, Fi and 
F 2 are the acoustic forces at the transducer faces, Z c is the acoustic impedance of the piezoelectric 
material„h 3 3 = e 33 /£ 33 is the piezoelectricity constant and C 0 is the capacitance value of the plate, and I 3 is 
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the electrical current. 0 = coejp/C D . The dephasage generates by the propagation with go = 2nf is the 
electrical pulsation where f is the frequency, e is the thickness of the piezoelectric plate, p is the material 
density, C D is the elasticity modulus with constant displacement field. 


Determination of the electrical impedance of the piezoelectric transducer 

The transducer is assumed here charged to the acoustic ports by homogeneous mediums. We 
determine the input electrical impedance of the transducer by the impedances of the acoustic loads Z\ and Z 2 
respectively acoustic impedance of the front medium and the back medium. At ports acoustic and given the 
sense of speed in the diagram, we have: 


F 1 = -Z x v x 

(3) 

F 2 = —Z 2 v x 

(4) 


We report these relations in the first two equations of the matrix (2) which makes it possible to 
determine the electrical impedance j- 


V 3 _ 1 1 _j_ C 0 h§ 3 -2Z c [l-cos 9\+j(Z 1 + £ 2 ) sin 9 

I 3 jC 0 co a) {Zq+Z 1 Z 2 )s\nO-jZ c {Z 1 +Z 2 )cos 6 


(5) 


The expression of electrical impedance is generally found in the literature in the following form: 
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( 6 ) 


K 2 q 1^2 

Where = 0 33 For longitudinal waves, the parameter K T is often defined as the piezoelectric 

9 CO 

coupling constant for a transversely clamped material; for it is the effective piezoelectric constant used when 
there is no motion transverse to the electric field. 


3. ELECTRO-ACOUSTIC MODEL OF PIEZOELECTRIC ELEMENT 

For reasons of simplicity, we chose to use Redwood's equivalent electrical circuit. The equivalent 
circuit of a thickness-mode piezoelectric transducer can be represented by the Redwood model as shown in 
Figure 2. 


Redwood model 



The model is divided in two parts. The first one is the electrical port which includes the capacitors 
C 0 and — C 0 that represent the capacitance motional effect. This electrical port is connected to a resistance R 
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and a voltage source noted V 3 .The second part is composed by the two acoustic ports, T is an ideal electro¬ 
acoustic transformer with a ratioh 33 C 0 . Piezoceramic layer is assimilated to a propagation line characterized 
by its characteristic impedance Z 0 (Z 0 = Z C .A), A is the area and the propagation timeTd (Td = e/v D ), with 
v D is the acoustic velocityand e is the thickness. 

One branch of the piezoceramic layer is in contact with the back medium Zback and the other is in 
contact with the propagation medium Zf ron t. 


4. VHDL-AMS BEHAVIORAL MODEL OF TRANDUCER 

Before discussing the modeling of transducers in VHDL-AMS, we quickly present this 
programming language. The VHDL-AMS is an IEEE standard [15]. It was developed as an extension of 
VHDL to enable the modeling and simulation of circuits and analog and digital - analog mixed systems. 

The VHDL-AMS implementation of the previous model as shown in Figure 2 is divided in two parts. First 
is the declaration of the entity which is composed of the physical characteristics of the transducer and the 
different terminals used in connection. Each TERMINAL depends of the physical nature of the relation to be 
implemented to describe the element. Using the statements ELECTRICAL in the electrical domain and 
KINEMATICV in the acoustic domain. The second part of the model is the architecture which establishes 
the physic laws related to the mathematical relation between each terminal [16]. 

In this approach, electrical and acoustic elements are described. Electrical element is defined with 
voltage and current quantities and acoustic element with force and velocity quantities [17, 18]. The 
VHDL-AMS implementation of Redwood's equivalent model as shown in Figure 2 is given by: 

ENTITY Redwood IS 
GENERIC (Co, kt, Z 0 , Td: real); 

PORT (TERMINAL p, m: electrical; TERMINAL tl 1, ml 1, t22, m22: kinematicv); 

END ENTITY Redwood; 

ARCHITECTURE structure OF Redwood IS 
TERMINAL pi: electrical; 

TERMINAL tl, tlx, t2x : kinematic v; 

QUANTITY vi across ii through p TO m; 

QUANTITY V 2 across i 2 through p TO pi; 

QUANTITY vte across ite through pi TO m; 

QUANTITY pti across uti through tl TO kinematicvground; 

QUANTITY plxr across ulxr through tlx TO tl; 

QUANTITY p2xr across u2xr through t2x TO tl; 

QUANTITY plx across u2x through tlx TO tl 1; 

QUANTITY p2x across ulx through t2x TO t22; 

QUANTITY pi 1 across tl 1 TO tl; 

QUANTITY p22 across t22 TO tl; 

BEGIN 
il==Co * vfdot; 
i2==-C 0 * v 2 'dot; 
pti==kt * vte; 

Uti==-ite/kt; 

p 1 xi==p22’DELAYED (Td) - plx; 
p2xr==p 11 ’DELAYED (Td) - p2x; 
p 1 x=(u 1 x+u2x’DELAYED (Td))*Z 0 /2.0; 
p2x=(u2x+u 1 x’DELAYED (Td))*Z 0 /2.0; 

END ARCHITECTURE structure; 


5. SIMULATION RESULTS 

We presentons in this part the results of simulation with VHDL-AMS of the input electrical 
impedance of the transducer. The studied transducer is built with PZT ceramic of PI-88 type, produced by 
Quartz and Silice society, vibrating in thickness mode at a frequency of 2.25 MHz, 1 mm thick and area 
A=132.73 mm 2 . These characteristics are recalled in Table 1. 

c 0 = £33 A /e ; v D = yfeiJp ; K t = h 33 j4JC° 3 ; Z c = pv D 
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Table 1. Transducer acoustic characteristics [19] 


Parametrs 

Definition 

Value 

Unit 

P 

Density 

7700 

Kg/m 3 

v D 

Acoustic velocity 

4530 

m/s 

Zo 

Acoustic impedance 

34.9 

Mrayls 

C 0 

Capacitor of the ceramic 

759 

PF 

K t 

coupling factor 

0.49 

- 

C33 

Elastic constant 

15.8xl0 10 

N/m 2 

fc 33 

dielectric constant 

870s 0 

F/m 

K33 

piezoelectric constant 

1.49xl0 9 

- 

tan 8 e 

Dielectric loss factor 

0.02 

- 


The frequential transducers response study is essential to predict the sensitivity of the system for the 
various analyzed mediums. Table 2 gives acoustic characteristics of simulated mediums. The amplitude of 
the transducer is fixed at lVolt with a frequency of 2.25MHz. 


Table 2, Acoustic characteristic of different mediums [20] 


Mediums 

Acoustic impedance (Mrayls) 

Acoustic velocity (kW) 

Water 

1.49 

1469 

Ethanol 

0.91 

1158 

Paraffin oil 

1.86 

1420 


5.1. Study of the electrical impedance of the Redwood model 

The complex impedance of the Redwood model is of the form: Z T (jw) = R + jX = |Z T |e J0 . We 
denote fi and f 2 the frequencies for which the impedance module reaches respectively its minimum and 
maximum values, the frequencies of resonance fr and antiresonance fa are the values for which the global 
reactance X of the transducer is zero (0=0). The simulation code used to obtain the transducer imput 
electrical impedance is given by: 

ENTITY ImpedanceSimulation IS 
END Impedance Simulation; 

ARCHITECTURE struct OF Impedance Simulation IS 
TERMINAL nl, n2: ELECTRICAL; 

TERMINAL n3, n4: KINEMATIC V; 

CONSTANT e: real:=1.0e-3 
CONSTANT A: real:=132.73e-3 
CONSTANT v D \ real:=4530.0; 

CONSTANT fo: real:=2.25e6 
CONSTANT Z c : real:=34.9e6; 

CONSTANT kt: real:=0.49; 

CONSTANT epsi 0 : real:=8.8542e-12; 

CONSTANT epsi 33 : real:=650.0; 

CONSTANT r 0 : real:=3300.0; 

CONSTANT h: real:=kt*y D *sqrt(ro/(epsio*epsi 33 )); 

CONSTANT Co: real:=A*epsi 33 /e; 

QUANTITY vac: real spectrum 1.0, 0.0; 

QUANTITY vinput across ie through nl to electricalground; 

BEGIN 

Vinput==vac; 

R: ENTITY resistanc (bhv) GENERIC MAP (50.0) PORT MAP (nl, n2); 

Tl: entity Redwood (bhv) generic map (Co, K, A*Z C , e/v D ) PORT MAP (n3, kinematicvground, n4, 
kinematicvground, n2, electrical ground); 

Rfront: ENTITY resistanc (bhv) GENERIC MAP (1.5e6*A) PORT MAP (n4, ground); 

Rback: ENTITY resistanc (bhv) GENERIC MAP (445.0 A*) PORT MAP (n3, electrical ground); 

END ARCHITECTURE struct; 

Figure 3 represents different characteristic curves that can be extracted from the impedance Z T (Jw ). 
In practice the frequencies Land fr on the one hand and the frequencies fa and f 2 on the other hand are 
very close. 
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Figure 3. Characteristic curves of the impedance of the Redwood model and Impedance phase 


5.2. Influence of the propagation medium on the impedance modulus 

The vibration modes of the piezoelectric transducers are strongly influenced by the medium in 
which they are immersed or by the structure of which they are integral. In a homogeneous liquid or gaseous 
fluid, the corresponding acoustic load modifies the vibration velocity of the transducer and dampens its 
resonance [21], Retrodiffused echoes to the transducer generate a reception current. In the case of a thin 
piezoelectric blade bonded to a solid structure under stress, the influence of the latter is preponderant. These 
considerations show that the electrical impedance of the transducer is very sensitive to its environment [22]. 
The transducer is loaded on the front face by three different mediums: water, ethanol and paraffin oil.Their 
acoustic characteristics are shown in Table 2. The simulation result is shown in the following figure. Figure 4 
clearly shows the influence of the medium on the electrical impedance modulus and precisely on Zmax. We 
note a concordance between the result obtained by VHDF-AMS and the experimental result presented in 
[23], so a good modeling by VHDF-AMS. 




Frequency (MHz) 


Figure 4. Transducer impedance modulus and the phase for different propagation mediums 


5.3. Influence of coupling factor on the impedance modulus 

From the formula (6) it is found that the electrical input impedance of the transducer is directly 
related to the coupling factor Kt. The following figure shows the variation of the impedance modulus for 
different coupling factors. The transducer is loaded on the front face by water as propagation medium. 

As shown in Figure 5, the difference between the resonant and antiresonance frequencies is strongly 
influenced by Kt. The influence of Kt is important on Zmax but negligible on Zmin. 
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Figure 5. Transducer impedance modulus and phase for different values of Kt 


5.4. Influence of intrinsic dielectric losses on the impedance modulus 

It should be noted that the form (6) of the electrical impedance, does not show the dielectric losses 
in the ceramic. These can be taken into account when considering a transmission line with losses. 

Electrical losses is well known to be the resistance leak of the ceramic capacitance considered as a 
resistance which depends on the frequency described by relation (7) [24, 25]. 


e coC 0 tan(5 e ) 

With kt the coupling factor in thickness mode, C 0 is the capacitance of the ceramic, tan(5 e ) is the 
losses factor and 0 ) the pulsation. 

Figure 6 shows a growth with R e of the maximum value Z Tmax of the electrical impedance at the 
antiresonance frequency, like we can notice a reduction of the phase (the inductive effect decreases in the 
resonance zone when R e increases). The antiresonance and resonance frequencies as well as the minimum 
value Z Tmin are independent of R e . 



Figure 6. Variations of the impedance modulus and the phase as a function of the dielectric losses 


6. CONCLUSION 

In this paper, a new approach to the ultrasonic transducer modeling system is presented. Our study 
begins with the presentation of an electric and electroacoustic model of the transducer. From this model we 
determine the analytical expression of the input electrical impedance. We modeled and simulated the input 
electrical impedance of the piezoelectric transducer with VHDL-AMS, as well as the influence of the various 
parameters such as the coupling factor, the propagation medium and the electrical losses on the modulus of 
the latter. We have shown that the frequency analysis of the electrical impedance can precisely locate the 
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resonance zone to control and stabilize the operating frequency of the system. In addition, the use of the 
VHDL-AMS language has the advantage of combining multi-physical domains, and indicates that the 
simulation of an ultrasonic detection device, comprising both electronic components and (electromechanical) 
transducers, is possible with VHDL-AMS. Finally, this study makes it possible to obtain a good 
representation of the state, the behavior and the performances in real time of the transducer. 
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